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H-Bases  I:  The  Foundation 


H.  Michael  Moller  and  Thomas  Sauer 


Abstract.  The  H-basis  concept  allows  an  investigation  of  multivariate 
polynomial  spaces  degree  by  degree.  In  this  paper,  we  mention  its  con- 
nection to  the  Grobner  basis  concept,  characterize  H-bases,  show  how  to 
construct  them,  and  present  a procedure  for  simplifying  polynomials  to 
their  normal  forms.  Applications  will  be  given  in  [8] . 


§1.  Introduction 

We  consider  II,  the  ring  of  polynomials  in  xi, . . . , xn  with  coefficients  from  an 
infinite  field  IK,  i.e.  II  = K[xi, . . . , xn],  and  the  subsets  II ^ of  all  polynomials 
of  degree  at  most  d.  In  many  applications,  one  is  interested  in  getting  a basis 
or  a generating  set  for  the  linear  vector  space  I Pi  lid,  where  I C II  is  an 
ideal.  Having  an  H-basis  {/i, . . . ,/3}  for  I,  then  the  set  of  all  p,  ■ /,  with 
Pi  € n d-deg(f i)>  i = generates  / fl Ify  as  a linear  vector  space.  Thus 

the  H-basis  concept  is  a tool  for  transforming  a non-linear  problem  in  n into 
a problem  in  one  (or  in  a series  of)  finite  dimensional  linear  space(s)  Ify. 

H-bases  were  introduced  first  by  Macaulay  [4].  His  original  motivation 
was  the  transformation  of  systems  of  polynomial  equations  into  simpler  ones. 
The  power  of  this  concept  was  not  really  understood,  presumably  because 
of  the  lack  of  facilities  for  symbolic  computations.  When  Computer  Algebra 
Systems  came  up,  Grobner  bases  (G-bases  for  short)  were  used  instead  of 
H-bases.  These  bases,  originally  invented  by  Buchberger  [2]  for  computing 
multiplication  tables  for  factor  rings,  are  now  also  applied  for  simplifying 
some  problems  in  Numerical  Analysis,  see  [5]. 

The  G-bases  give  generating  systems  not  to  I fl  but  to  I fl  Ti,  where 
T{  C n is  a linear  vector  space  of  dimension  i,  and  Ti  C T+i  for  all  i and 
n = U,>o^.  This  finer  decomposition  has  some  drawbacks.  For  instance  if 
an  ideal  is  invariant  under  an  affine  symmetry  group,  its  G-bases  are  typically 
not  invariant.  Since  the  spaces  II  d are  invariant  under  affine  symmetry  groups, 
H-bases  do  not  destroy  such  symmetries. 

Many  of  the  problems  in  applications  which  can  be  solved  by  Grobner 
techniques  can  also  be  treated  successfully  with  H-bases.  In  [7]  we  gave  an 
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overview  of  such  problems.  In  the  present  paper,  we  describe  briefly  the 
underlying  concept  of  grading  rings,  which  leads  to  G-  and  H-bases,  and 
present  some  properties  characterizing  H-bases.  In  contrast  to  [7],  where  we 
only  gave  a class  of  examples  of  H-bases,  we  present  here  the  construction  of 
an  H-basis  for  zero-dimensional  ideals  I.  A useful  tool  for  our  procedure  is 
the  so  called  normal  form  mapping  NF,  presented  in  Section  4,  which  projects 
n orthogonally  to  the  ideal  I provided  an  H-basis  of  / is  given.  In  [8]  we  show 
how  these  normal  forms  can  be  applied  in  numerical  applications. 

§2.  H-bases  and  G-bases 

In  ring  theory,  rings  can  be  graded  by  an  ordered  monoid,  i.e.  by  an  abelian 
semigroup  F with  addition  + and  total  ordering  -<  satisfying 

7i  -<  72  =>  7o  + 7i  -<  7o  + 72,  V70, 7i> 72  € T. 

There  are  two  major  examples  for  grading  n by  an  ordered  monoid  F : 

n = © <),  n£>n£>  c n<^72  v7i , 72  € r. 

7«r 

The  first  one  is  the  H-grading  with  T :=  1N0, 

n(,r)  :=  {p  € n | p homogeneous  of  order  7}. 

The  ordering  of  T = 1N0  is  the  natural  one.  The  second  example  for  gradings 
is  the  G-grading,  where  T :=  IN"  and 

n(7i,...,7„)  :=  W'  ' " Xn"  I C G K)' 

T = IN"  is  ordered  by  an  admissible  term  ordering, 

Odt,  i<j=>i  + k-<j  + k. 

Since  the  decomposition  of  n into  the  sets  II  ^ is  a direct  sum,  every  f e H 
has  a unique  representation  f = YL  f-r  The  maximal  7 with  /7  0 is  called 

the  maximal  part  of  / / 0,  for  short.  It  is  also  called  the  maximal  form 

in  the  H-case,  or  leading  monomial  in  the  G-case.  In  the  G-case,  = 

lc(f)lt(f ),  where  lc(f)  € K is  the  leading  coefficient  and  = x'f  ■ ■ • x]? 

the  leading  term.  The  maximal  form  of  / 0 is  also  denoted  by  M //(/). 

Definition  1.  {pi,  ■ ■ ■ ,pm}  C I is  called  a basis  of  an  ideal  I C n,  briefly 
I = 'TV  pel 

m 

351, •••,5m  e n : P = Y^gkPk- 

k= 1 
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It  is  also  a G-basis  or  H-basis  resp.  if  <?i, . . . , gm  satisfy  in  addition 

max  lt(gk)lt(pk)  = lt(p)  (G  — basis), 
k= 1 


or  max  deg(gkpk)  = deg(p)  (H  - basis). 

k—1 

Theorem  1.  Let  I = (pi, . . . ,pm).  Then  {pi, . . . ,pm)  is  an  H-basis  (G- 
ba sis  resp.)  if  and  only  if  the  least  ideal  containing  all  Mff(f),  Of - f £ I 
(or  all  lt(f),  0 / / 6 / resp.)  is  generated  by  MH(pi), . . .,MH(pm)  (or  by 
lt(pi),  ■ ■ ■ ,lt(pm)  resp.). 

This  theorem,  which  holds  mutatis  mutandis  for  arbitrary  graded  rings, 
is  proved  for  instance  in  [6].  An  immediate  consequence  of  it  is  that  every 
ideal  I (0)  has  an  H-  and  a G-basis. 

G-bases  are  now  a standard  tool  in  Computer  Algebra.  They  are  covered 
by  nearly  all  textbooks,  and  are  contained  in  almost  all  Computer  Algebra 
Systems.  The  grading  by  one-dimensional  linear  spaces  II^P  often  reduces  the 
computation  to  solving  a series  of  one-dimensional  problems.  On  the  other 
hand,  the  construction  of  G-bases  is  often  difficult  or  even  impossible  because 
of  the  high  complexity  of  Buchberger’s  algorithm  for  computing  G-bases.  In 
addition,  in  many  applications  the  G-bases  allows  only  little  insight  into  the 
structure  of  a solution  by  the  artificial  ordering  term  by  term. 

§3.  Characterization  of  H-bases  and  Normal  Forms 

Macaulay  introduced  H-bases  using  homogenizations  and  dehomogenizations 
of  polynomials.  The  name  H-basis  originates  from  the  first  letter  of  homoge- 
nization. 

Definition  2.  Let  f € K[xj, . . . , xn]  have  degree  d, 

d 

f = fi  homogeneous  of  degree  i,  fd  ff  0. 

i=0 

Then  introducing  a new  variable  xo,  the  homogenization  of  f is  a homogeneous 
degree  d polynomial  in  K[xo,  x\, ... , xn], 

m :=X>o-\A- 

i= o 

A homogeneous  F € K[a;o,  ®i, . . . , xn\  can  be  dehomogenized  to  an  f £ U by 

Xq  = 1. 

For  more  details  on  homogenizations  and  their  connection  to  projective 
coordinates,  we  refer  to  [3J. 
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Theorem  2.  (Macaulay  [4]}.  Let  I = (hi,...,  hs).  Then  the  following  state- 
ments are  equivalent 

1)  The  least  ideal  containing  all  $(/i),  0 ^ h £ I,  is  (<f>(hi), . . . , $(hs)). 

2)  X0F  e , *(/».)>  =>  F e ($(/n), . . . , $(/»,))• 

3)  {/ii, . . . , /is}  is  an  H-basis  of  I. 

The  power  of  the  G-basis  concept  is  mainly  based  on  the  possibility  of 
reducing  a polynomial  to  a simpler  one  by  subtracting  suitable  multiples  of 
elements  of  the  G-basis.  A consequent  application  of  this  reduction  strategy 
gives  the  so  called  normal  form,  in  a sense  the  simplest  polynomial  obtainable 
by  the  reductions.  We  translated  this  technique  to  H-bases  in  [7],  and  give 
here  for  consistency  a short  resume  of  the  main  results. 

Definition  3.  We  denote  by  the  space  of  all  homogeneous  degree  d 

polynomials  for  d £ 1N0  and  II ^ :=  {0}  for  d < 0.  Let  hi,. . . ,hs  £ II.  Then 
we  define  a finite  dimensional  linear  subspace  of  II by 

S 

vd(hi,. ..,hs)  :=  1 9i  € n^de9(hi)}. 

i=l 


Analogously  for  an  ideal  1 C II, 

Vd(I)  :=  {M„(p)  | p £ /,  deg(p)  = d)  U {0}. 

We  introduce  an  inner  product  ( . , .)  in  IT,  for  instance,  by  the  inner 
product  of  the  (weighted)  coefficient  vectors,  or  by  a strictly  positive  linear 
functional  J and  (/,  g)  :=  J(f  ■ g)  if  K C R or  :=  J(fg ) if  K = C.  Then  we 
can  define  orthogonal  complements  Wd(hi, . . . , hs)  and  Wd(I)  in  . Hence 

Vd(hi,. . . , hs)  ® Wd{hu. . . , hs)  = U{dH)  and  Vd(I)  ® Wd(I)  = U{dH) . 

Let  us  consider  a polynomial  / of  degree  d.  Then 

€ Vd(hi, . . . , hs)  © Wd(hu. . . , hs). 

Let  wd  denote  its  natural  projection  on  Wd(hi, . . . , h$).  This  homogeneous 
polynomial  can  be  computed  by  solving  a finite  linear  system  of  equations  be- 
cause  nd  ; has  a finite  dimension.  Hence  there  are  homogeneous  polynomials 
gi, ...  ,gs  such  that 

S 

f = wd  + '^2,gihi  + /1,  gi  € /1  £ nd_i. 

1=1 


We  say  / reduces  to  wd+f\  modulo  {hi, . . . , hs}  and  call  f\  then  the  remainder 
of  /. 
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In  the  reduction  modulo  {hi, . . . , hs}  the  degree  of  the  remainder  fi  is  less 
than  deg(f).  Hence  this  reduction  can  be  applied  recursively  reducing  /*_ 1 
constructively  to  wj+i-i  + fi  modulo  {hi, . . . , hs}  for  i = 1, . . . , d+ 1 starting 
with  / = /o  £ n<j  and  terminating  with  fd+i  = 0,  because  the  constant  f,i 
is  either  in  Vo(hi, . . . , hs)  or  in  Wo(hi,  • • • , hs).  Combining  these  reductions 
modulo  {hi, . . . , hs},  one  obtains  for  / 

/ = ^ Wi  + Y'Yjijhj,  gij  £ n|_]e3(ft.r 

t=0  i=0  j— 1 

Then  J2f=owi  uniquely  determined  by  /,  by  {hi, . . . ,hs},  and  by  the  un- 
derlying inner  product. 

Definition  4.  Let  hi,...,hs  £ II.  We  say  f £ Hd  reduces  fully  modulo 
{hi,...,hg}  to  Ya= owi  every  w,  £ Wi(hi, . . . ,hs)  is  constructed  as  de- 
scribed above.  Ylj=owj  JS  called  the  normal  form  of  / modulo  {hi, . . . , hs}, 
for  short 

d 

NF  (/,  {hi,...,hs})  :=J2WL 
3=0 

If  {hi, . . . , hs } is  not  an  H-basis  of  I :=  (h i, . . . , hs),  then  is  not 

necessarily  contained  in  Vie9(/)(hi, . . . , hs),  although  / £ I.  This  means,  that 
eventually  the  first  homogeneous  polynomial  Wd  is  not  0 if  / £ I.  Hence  at 
most  if  {hi, . . . , hs}  is  an  H-basis,  then  NF (/,  {hi, . . . , hs})  = 0.  In  fact,  as 
quoted  in  [7]  but  shown  already  in  [9],  {hi,. . . ,h5}  is  an  H-basis  if  and  only 
if  NF(/,{hi,...,hs})  = 0 for  every  f £ (hi, . . . , hs). 

Another  characterization  of  H-bases  given  in  [7]  is  as  follows. 

Theorem  3.  Let  I be  an  ideal  and  hi, . . . ,hs  £ I.  {hi, . . . , hs}  is  an  H-basis 
of  I if,  and  only  if,  for  all  d £ IN, 

Vd(I)  = Vd(hu...,hs). 


§4.  On  the  Construction  of  H-bases 

Macaulay  proposed  in  [4]  a procedure  for  computing  H-bases  of  ideals  given  by 
a basis.  However,  his  description  was  only  by  an  example.  He  claimed  “ This 
procedure  is  a general  one ” . But  he  needs  in  his  example  the  computation  of 
certain  modules  of  syzygies.  These  can  be  constructed  only  in  special  cases 
or  by  computing  first  a G-basis  and  then  applying  techniques  as  in  [1], 

On  the  other  hand,  if  an  admissible  term  ordering  -<  is  compatible  with 
degrees, 

degixj1  ■ ■ ■ xln)  < deg(x^  ■ ■ ■ ar£”)  =>  xj1  • ■ ■ x%n  ■<  x^1  ■ ■ ■ x%n , 

then  a G-basis  with  respect  to  -<  is  also  an  H-basis.  Hence  Buchberger’s 
algorithm  for  computing  G-bases  also  serves  for  computing  H-bases.  This 
seems  a more  direct  access  than  via  syzygies.  However,  if  one  wants  to  use 
H-bases  instead  of  G-bases,  this  way  is  still  a detour. 
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In  case  the  number  n of  variables  coincides  with  the  number  of  given 
polynomials,  then  there  is  an  easy  test  for  H-bases  as  proved  in  [7]. 

Theorem  4.  Let  h\ , . . . , h„  be  n polynomials  such  that  their  maximal  forms 
Mu{hi), . . . , MH{hn)  have  only  the  point  (0, . . . , 0)  as  common  zero.  Then 
{hi, . . . , hn } is  an  H-basis. 

For  an  arbitrary  zero-dimensional  ideal,  i.e.  for  an  ideal  I such  that  the 
polynomials  in  7 ( equivalently:  the  polynomials  in  an  arbitrary  basis  of  7) 
have  only  a finite  number  of  common  zeros  in  K , K the  algebraic  closure  of 
K,  see  [3],  we  present  here  a procedure  which  computes  an  H-basis  from  a 
given  basis. 

Procedure  for  computing  H-bases. 

In  : TL0,  a finite  polynomial  set  generating  a zero-dimensional  ideal  7. 
Out  : 77,  an  H-basis  of  7. 

Start:  ft  :=  ft0,  d = 0. 

Loop:  Check  the  finite  dimensional  linear  vector  space 

Vd(H)  :=  {Y^auMnih)  \ ghMH(h)  e n<ff)} 

hen 

for  linear  dependencies.  If  YLhenQ^iiih)  = 0,  then  compute  p := 
NF(£h6W  <?/»h,77)).  If  p ^ 0,  then  enlarge  77  by  p,  and  modify  con- 
sequently Vo(H), . . . , Lower  d to  the  least  k where  14(77)  is 

changed  and  go  to  Loop.  If  for  no  linear  dependency  such  p is  nonzero, 
then  then  enlarge  d by  1.  If  now  14(77)  = holds  true,  then  return 
ft  otherwise  go  to  Loop. 

This  informal  description  can  be  extended  easily  to  a correct  algorithm.  One 
has  to  observe  that  the  checking  of  Vd(ft)  for  linear  dependencies  needs  a basis 
of  the  nullspace 

S 

{(St  • ■ • ift)  £ x ...  x | y^ygjMH{hj)  = 0}, 

i=l 

where  ft  = {hi, . . . , hs}.  If  for  every  basis  element  {gi,...,gs)  the  normal 
form  of  gihi  is  0,  then  it  holds  for  every  element  of  the  nullspace,  i.e. 
for  every  dependency.  As  a byproduct  of  the  basis  computation  one  obtains 
dimVd{ft ).  Then  the  test  Vd{ft)  = reduces  to  a comparison  of  the 

dimensions  because  of  Vd{ft)  C . 

For  proving  correctness  and  termination,  we  consider  first  / :=  Ylhen  dh h 
with  ghMjj{h)  € U{dH)  for  all  h € ft.  If  NF(/,  ft)  = 0,  then  especially  Mjj{f)  € 
14(77)  for  a k < d,  and  hence 

MH(f)  € (MH(hi), MH(hs)),  where  ft  = {hi, . . . , hs}. 

In  case  p :=  NF(/,  77)  7^  0 either  M//(/)  7^  Mh{p)  holds,  i.e.  again 
Mn{f)  G ( MH{h\ ), . . . , MH(hs)), 
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or  Mn(f)  = Mh{p)  holds,  i.e. 

Mji(f)  G ... , MH{hs),  Mh{p))- 

Therefore,  if  in  the  procedure  d is  increased  by  1 (and  7i  is  updated),  then  for 
every  0 i=-  J2hen  9h^  e 7,  deg(gh)  + deg(h)  < d the  relation 

Mh(Y,  9hh ) € (Mnih), ...,  MH{hs )) 
hen 

holds  where  again  77  =:  {hi, . . . , hs}.  This  is  our  inductive  hypothesis. 

The  ideal  I has  an  H-basis,  say  {ipi, . . . , <pm}.  77  is  a basis  of  I.  Hence 
every  ipi  has  a representation  tp%  = 2j=i  9ijhj,  9ij  € n.  If  the  inductive 
hypothesis  holds  for  d,  then  one  obtains  at  least  for  d > M :=ma Xjjdeg(gijhj) 
that 

Mnivi)  6 (Mn(hi), ...,  MH(hs)),  i = 1, . . . , m. 

Hence  for  those  d Vd(ipi, . . . ,(pm)  C Vd(U).  But  Vd(ipi, . . . , <pm)  = Vd(I), 
since  {yi,. . . , Pm}  is  an  H-basis  of  I.  Therefore  Vd(H)  = Vd(I)  for  d > M. 
By  the  inductive  hypothesis,  also  Vk(H)  = 14(7)  holds  for  k < M.  Hence  H 
is  an  H-basis  of  I if  we  arrived  at  a d > M in  the  procedure. 

The  ideal  7 has  dimension  0.  Then  there  is  a D such  that  Vd(I)  — IlJ/^ 
for  all  d > D,  see  for  instance  [3,  Ch  9.4, Prop. 6]  and  [3,  Ch  5.3,Thm.6]. 
Hence  for  d > max{D,  M}  one  has  14(77)  = n ^ . Thus  in  the  course  of  the 
procedure,  one  arrives  once,  not  knowing  M,  at  a do  with  Vdo{TL)  = 

Then  also  14(77)  = II for  all  k > do.  Therefore,  for  every  polynomial  / 6 7 
with 

5 

f = y Zsihj,  g{  € n{kJdeg{hi) 

i= 1 

the  assumption  ^ (M/f(hi), . . . , Mh(/is))  leads  to  deg(Mu(f))  < do- 

But  then  the  inductive  hypothesis  gives  a contradiction.  Therefore  the  proce- 
dure gives  no  new  p ^ 0 enlarging  the  set  77.  This  ensures  termination. 

An  implementation  of  an  algorithm  based  on  this  procedure  and  a com- 
plexity analysis  is  still  a work  under  progress. 
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